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V. V. Skokov 1 **, V. D. Toneev** 

* Bogoliubov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, 141980, Dubna, Russia 
' Gesellschaft fur Schwerionenforschung, 64291 Darmstadt, Germany 

A hybrid model is put forward for describing relativistic heavy ion collisions. The early inter- 
action stage responsible for entropy creation is calculated within the transport Quark-Gluon 
String Model resulting in an initial state. The passage to subsequent isoentropic expansion 
proceeds with exact account for all conservation laws. Relativistic 3D hydrodynamics with 
the mixed phase equation of state is applied to this expansion stage. Essential differences 
in evolution of a fireball, described within this model and one for the Bjorken regime, are 
noted. The in-medium modified e + e _ spectra are studied and confronted with the recent 
CERES/NA45 data for 8% central Pb+Au collisions at the bombarding energy 158 AGeV. 

PACS: 25.75.-q24.10.Nz24.10.Pa 

Electromagnetic probes and, in particular, di-electrons play an exceptional role providing al- 
most undisturbed information on a state of highly compressed and hot nuclear matter (fireball) 
formed in relativistic heavy ion collisions. Generally, dilepton yield depends on both global prop- 
erties of matter constituents (hadrons and/or quarks, gluons) defined by the equation of state and 
also on individual constituent properties related to their in-medium modification. Such a mod- 
ification has been observed by the CERES Collaboration in the analysis of the e + e~ invariant 
mass spectra from central Pb+Au collisions [1]. The measured dilepton excess in the range of in- 
variant dilepton masses 0.2 < M < 0.7 GeV may be interpreted in terms of a strong in-medium 
p-meson modification (see review articles [2, 3]). Many papers are devoted to the analysis of 
these CERES data but mainly two scenarios of hadron modification are available on the market: 
Consideration based on the Brown-Rho scaling hypothesis [4] assuming a dropping p mass and 
that for a strong broadening as found in the many -body approach by Rapp and Wambach [3, 5]. 
Unfortunately, a very simplified dynamical treatment is usually applied which obscures the origin 
of the model agreement or disagreement with experiment. 

In this paper, we make emphasis on nuclear interaction dynamics to clarify which states in 
the main contribute to the observed dilepton yield as well as which are the model parameters 
used. We confront our model calculations with the recent CERES/NA45 results for 8% central 
Pb+Au collisions at the bombarding energy 158 AGeV [6], 

We treat dynamics of heavy ion collisions in terms of a hybrid model where the initial inter- 
action stage is described by the transport Quark Gluon String Model (QGSM) [7]. 
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Fig. 1. Average entropy S per average baryon charge Qb for participants versus relative distance between 
the centers of two colliding Pb nuclei at Ei a t =158 AGeV. 



In FigQ the temporal dependence of the ratio of entropy S to the total participant baryon 
charge Qb is shown for Pb+Pb collisions at the impact parameter b — 2.5 fm and bombarding 
energy 158 ^4GeV evaluated within QGSM. Being calculated on a large 3D grid, this ratio is less 
sensitive to particle fluctuation as compared to entropy itself. Small values of Qb at the very 
beginning of interaction result in large values of the S/Qb ratio. It is clearly seen that this ratio 
is practically kept constant for R > 1.5 fm, i.e., a little bit later on the moment when two nuclei 
completely overlap (R = 0). The stage for R> 1.5 fm may be treated as isoentropic expansion. 

So the subsequent stage starts from R w 1.5 fm, which corresponds to the kinetic evolu- 
tion time moment tki n = 1.8 fm/c and evolves further as a isoentropic expansion. This stage 
is described within the relativistic 3D hydrodynamics. A full set of hydrodynamic equations is 
solved numerically by the FCT-SHASTA method [8] with the mixed phase Equation of State 
(EoS) [9]. This thermodynamic ally consistent EoS uses the modified Zimanyi mean-field in- 
teraction for hadrons which reproduces properties of the saturated nuclear matter at the normal 
density tlb = no and is in good agreement with the Danielewicz constraints to EoS [10] (see 
Fig0. This mixed phase EoS also includes the interaction between hadron and quark-gluon 
phases, which results in a cross-over deconfinement phase transition. In addition to the model 
prescription in Ref. [9], the hard thermal loop term was self-consistently added to the interac- 
tion of quarks and gluons to get the correct asymptotics at T >> T c and to provide reasonable 
agreement between the model results and lattice QCD calculations at finite temperature T and 
chemical potential ub [11]. 

To proceed from kinetics to hydrodynamics keeping exactly the conservation laws, we eval- 
uate locally on the 3D grid the energy-momentum tensor, baryon density and velocity profile at 
the time moment tun = 1.8 fm/c and then treat these quantities as an initial (t=0) state of a 
fireball for its subsequent 3D hydrodynamic evolution. The temporal dependence of the average 
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Fig. 2. Baryon density dependence of pressure at the zero temperature. The mixed phase model results are 
plotted by the solid line, the dashed one corresponds to the ideal gas EoS. The shaded area is the constraint 
obtained by Danielewicz et al. [10]. 



energy density is presented in Figs|3]for the expansion stage. The average quantity is defined as 

£ -( t ) S£ -=/^ e( «,x)»(«,x)//^»(M), 

where w(t, x) is its time-space distribution. 

As is seen, the initial energy density is quite high, e(t w 0) ~ 5.5 GeV/ fm 3 . The appropri- 
ate compression ratio ub (0)/ no reaches about 4.2. Note that these values characterize the initial 
locally equilibrium state of the fireball and are noticeably lower than the peak-values estimated 
in nonequilibrium models with averaging over some limited spacial volume [13]. Average en- 
ergy density falls down in time and its behavior at the SPS energy remarkably differs from that 
for the Lorentz invariant Bjorken expansion regime (see the dashed line in FigQ. The scaled 
Bjorken solution of ID hydrodynamics [12] is not applicable to the first moments of hydrody- 
namic expansion and more corresponds to the asymptotic expansion regime. However, if the 
same initial energy density is assumed for the Bjorken expansion, one can see that the fireball 
life time becomes very short. In addition, this simple solution does not reproduce some flattening 
of the time dependence of the energy density observed for t > 8 fm/c. The origin of this effect 
is demonstrated in the right side of Fig|5] 

In our consideration, the hydrodynamic evolution of a cell ends at the freeze-out defined 
by the condition that the local energy density (including 6 neighbor cells) is below a certain 
value Si < Ef r = 0.3 GeV/fm 3 . This condition may be fulfilled locally at the periphery of 
the expanding system from the very beginning of expansion and works continuously during the 
whole evolution. The frozen cell will not participate in a further dynamical process but contribute 
to the "hadron decay cocktail". The system volume V presented in Fig|3]is defined as the sum 
over all cells with nonvanishing energy density £j besides frozen cells. As is seen, the system 
volume naturally starts to grow with expansion but later on the competition between expansion 
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Fig. 3. The average energy density (the left panel) and volume (the right panel) of the expanding fireball 
formed in central Pb+Pb collisions. The dotted line (in arbitrary normalization) shows the Bjorken solu- 
tion [12] for ID expansion with the ultrarelativistic ideal gas EoS and the dot-dashed line corresponds to 
the account for hadronic phase cells only (see the text). 



and freeze-out effect is getting stronger and so V drops for t > 8 fm/c. Certainly, this effect is 
absent in the Bjorken regime. According to the mixed phase model [9], in every state generally 
there is a homogeneous mixture of free quarks/gluons and bound quarks (hadrons) with quark 
densities p®^~ arks and p^ uar y., respectively. One may conditionally define the hadronic phase as 

P^uark > Pquarks- ^ s f°ll ows from Fig|3](the dot-dashed line), the admixture of the quark/gluon 
phase is negligible for the expansion time t > 4 fm/c which roughly corresponds to the average 
temperature of the system T(t w 4) « 180 MeV. 

The full set of the freeze-out cells or the time-space freeze-out surface, where hadrons are 
decoupled, defines all characteristics of observable hadrons. To transform these frozen cells into 
hadrons, we use the procedure described in [14]. The freeze-out parameter £f r is controlled by 
the pion yield at high colliding energies (but weakly sensitive at moderate energies, similarly 
to [14]). Having in mind that the main emphasis will be made on the irir channel, in Fig@]the 
rapidity and transverse mass spectra of observed pions are shown for central Pb+Pb collisions at 
the top SPS energy. The agreement between model calculations and data of the NA49 experiment 
is quite reasonable. 

To find observable invariant mass spectra of dileptons, one should integrate the emission 
rate d 8 N ee (T(x), /j,b{x), M,rj,qT)/dt d 3 x d 4 q over the whole time-space evolution J dtd 3 x 
and add the contribution from the freeze-out surface (hadron cocktail). If we are interested in the 
e + e~ invariant mass M distributions in some pseudo-rapidity window A?7 e ± (like in the CERES 
experiment) and treat the fireball evolution on average, then in this case 

d 2 N ee M f f [ 2 « [°° d 5 N ee (f(t),fL B (t),M, v ,g T ) 

dm = A^J dv J vm Jo d *L qrdqT **h Ac( 
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where the factor Acc = Acc(M, 77, qr) takes into account the experimental acceptance. 

To simplify more our task, we consider only the main emission channel irir — > . The 
dilepton emission rate is [16] 

d 5 N , s a 2 . . . — , 

= -C(M) -=-3 f B (qo,T) lmU em (q,T,fl b ), (2) 



dtd^q TT 3 q 

where the Bose distribution function is defined as /s(?oj T) = [e q °l T — l) -1 , the 4-momentum 
transfer q 2 = M 2 = q 2 — q 2 and the lepton kinematic factor is 



with the lepton mass m; . 

We will consider electro-magnetic current correlation function H em (q, T, fib) in two different 
scenarios: (I) the finite temperature modification of pion gas properties [16]; (II) droping mass 
scenario [3]. 

Dilepton invariant mass spectra from 8% centrality Pb+Pb collisions at Ei a b = 158 AGeV 
are presented in Fig[5] Experimental acceptance is taken into account for both hadron cocktail 
and 7T7T annihilation channel. This channel is calculated on the absolute basis without any fitting 
parameters. As one can see, the 1st scenario slightly overestimates the dilepton yield near the 
p meson mass and underestimates it in the 0.3 < M < 0.6 range though these values are es- 
sentially above the hadron cocktail data. This overestimation is explained mainly by simplified 
dynamics of dilepton production, using average evolution trajectories and rough approximations 
for proceeding to hadronic phase. The depleted yield of low-mass dileptons results from neglect- 
ing "sobar" contribution (nucleonic loops) into the emission rate. 
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Fig. 5. Invariant mass distribution of dileptons from 8% central Pb+Pb collisions at the beam energy 
158 AGeV. Experimental points are from [6]. The solid and dashed curves are calculated using the finite 
temperature modification in pion gas and the p dropping mass scenario J4}, respectively. The dotted line 
indicates an estimate of "hadron decay cocktail". 



To simulate droping mass scenario, we modify the in-medium mass of p-meson as commonly 
employed in the literature : 

m* p = m p (l - 0.18 • n B {t)/n ) (l - [T(i)/T c ] 2 ) ' 3 (4) 

and simultaneously apply the same modification to the vector dominance coupling m 2 p /g — > 
m* 2 /g*. The first factor in eq.@ was estimated by Hatsuda and Lee [18] according to the QCD 
sum rules. The temperature dependence in eq.@ is motivated by the T-dependence of the quark 
condensate. 

The results for the dropping p mass scenario are also plotted in Fig[5] In this case, the low 
dilepton mass region 0.2 < M < 0.7 is also populated but the free p mass region is essentially 
underestimated. The shape of the ill-distributions in these two scenarios is quite different. The 
simplified dynamical assumptions mentioned above also influence the p mass scaling results. 
However, the large shift of the e + e~ invariant mass spectrum towards low M is mostly due to 
using the T-dependent factor in eq.@, which seems to be not very justified [19]. 

Finally, inclusion of the full dynamics of heavy ion collisions in evaluation of dilepton pro- 
duction is quite important. Comparison of evolution within the proposed dynamical model with 
the popular Bjorken regime shows essential differences which may manifest themselves in dilep- 
ton observables. Full dynamical calculations allow one to probe not only the shape of the e + e~ 
spectra from 7T7t annihilation but also to estimate the absolute contribution of this channel. Cer- 
tainly, the simplifying dynamical assumptions used should be overcome as well as other dilepton 
production channels should be added, in particular those coming from a quark phase. In this re- 
spect, it is of interest to note that the mixed phase EoS, in principle, opens new possible channels 
due to interaction between hadron and quark/gluon phases. The hadron decay cocktail should 
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also be calculated within the same model. This work is now in progress. 
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